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; 1 INTRODUCTION 



ABSTRACT 

We extend the Hubble diagram up to z = 5.6 using 63 gamma-ray bursts (GRBs) 
via peak energy-peak luminosity relation (so called Yonetoku relation), and obtain con- 
straints on cosmological parameters including dynamical dark energy parametrized by 
P/p = w(z) = wq + w a ■ z/(l + z). It is found that the current GRB data are consis- 
tent with the concordance model, (fi m = 0.28, Q\ = 0.72, wo — —l,w a — 0), within 
two sigma level. Although constraints from GRBs themselves are not so strong, they 
can improve the conventional constraints from SNela because GRBs have much higher 
redshifts. Further we estimate the constraints on the dark-energy parameters expected 
by future observations with GLAST (Gamma-ray Large Area Space Telescope) and 
Swift by Monte-Carlo simulation. Constraints would improve substantially with an- 
other 150 GRBs. 
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In our previous paper (|Kodama et aljfeoOct ). we calibrated 
the peak energy-peak luminosity relation of GRBs with 33 
nearby events (z < 1.62) whose luminosity dis tances were 
estimated from those of large amount of SNela ( Riess et alj 
120071 : IWood-Vasev et all 120071 ; iDavis et all 120071 1. This cal- 
ibrated Yonetoku relation, derived without assuming any 
cosmological models, can be used as a new cosmic distance 
ladder toward higher redshifts. Then we determined the lu- 
minosity distances of 30 GRBs in 1.8 < z < 5.6 using 
the calibrated relation and calculated the likelihood varying 
(fi m ,fi A ). Weobtained (O m ,fi A ) = (0.371^, 0.631^) for 
a flat universe, which is consistent with the concordance cos- 
mological model within one sigma level. Our logic to obtain 
a new distance ladder is similar to that for SNela, that is, 
we calibrate a new distance indicator (Yonetoku relation) 
at low redshifts (z < 1.62) using the well established indi- 
cators (SNela). Then we assume the new relation holds at 
high redshifts (z > 1.62), although more detailed analysis 
is needed for possib le selection bias an d evolution effects in 
the relation (|Oguri fc Takahashi|[2006l '). 
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Currently the number of GRBs with z > 1.62 is rela- 
tively small (~ 30) and the statistical and systematic errors 
of the Yonetoku relation are not so small compared with 
SNela. Even then, GRBs are still effective to probe dark 
energy, especially for dynamical dark energy whose energy 
density becomes large at high redshifts, beca use the mean 
redshift of GRBs is higher than that of SNela (|Amati et al.l 
20081: iLiang et al.ll2008l; iBasilakos fc Perivolaropoulos!l200a 
Schaefeij|2007i : iGhirlanda et al.ll2006l : iFirmani et al.l 120060 . 



There would be many ways to characterize the time vari- 
ation of dark ene rgy. Here we adopt a simple phenomeno - 
logical model as (|Chevallier fc Polarskil [ioOll : iLinded l2003h 
P/p = w(z) = Wo +w a ■ z/(l + z) = Wo + w a (l — a), where 
wo and w a are constants and a is the scale factor of the uni- 
verse. For this model, GRBs would give strong constraints 
on w a which represents the time dependence of dark energy. 

In this Letter, we present constraints on cosmological 
parameters such as wo and w a using both SNela with z < 1.8 
and GRBs with z < 5.6. In § 2, we briefly review the main 
results of the previous paper (|Kodama et al.ll2008h and the 
theoretical and observational basis of the analysis. In § 3- 
1, we assume cosmological constant (w — —1) and obtain 
constraints on (il m ,QA)- In § 3-2, we assume a flat uni- 
verse and non-dynamical (w a = 0) dark energy and obtain 
constraints on (wq, Q m ). In § 3-3, we fix Q m = 0.28 for sim- 
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plicity, and obtain the plausible values of Wo and w a - In 
§ 4, we discuss how these constraints will be improved in fu- 
ture observations of high redshift GRBs by such as GLAST. 
Throughout the paper, we fix the current Hubble parameter 
as H = 66 km s _1 Mpc _1 . 



2 FRIEDMANN UNIVERSE 



In our previous paper dKodama et al-lbOOST ). we calibrated 
the Yonetoku relation IjYonetoku et alj2004f ) using 33 events 
with the redshift z < 1.62 as 



-) = (1.31 ±0.67) x 10" 4 ( 



E, 



V \1.68±0.09 
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10 52 erg s _1 ' 

where L p and E p are the peak luminosity and the peak en- 
ergy of the spectrum of a certain GRB event in a comov- 
ing frame, respectively. Fig. 1 of the previous paper shows 
that the linear correlation coefficient of the above relation 
in logarithmic scale is 0.9478 and the chance probability is 
6.0 x 10~ 17 . However the data distribution has a larger de- 
viation around the best fit line compared with the expected 
Gaussian distribution. We estimated this systematic devia- 
tion in the normalization as 9.57 x 10~ 5 . 

Now if one assumes that the above relation holds even 
for z > 1.8, we can determine the luminosity distance cLl{z) 
from the observed peak flux (/ Pj0 bs) and the observed peak 
energy of the spectrum E Pi0 b B since L p — 4tvcIl(z) 2 / p ,obs and 
E p — (l+z)E Pt obs. We express unknown equation of state of 
dark energy as P = wx(z)p and Six as the present energy 
density of the dark energy divided by the critical density. 
Then in the Friedmann universe with Qt = Sl m + Six — 1, 
the luminosity distance is given by 
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Using the distance modulus (fio(zi)) from observation of 
SNela and GRBs, we define a likelihood function as 
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Figure 1. Hubble diagram extended to z = 5.6 from GRBs. The 
green points and red points are th e luminosity distance deter- 
mine d by SNela jtliess et al] 12007ft and GRBs jKodama et al] 
2008), respectively. The blue line is the luminosity distance 
of ACDM model with (f2 m ,n A )=(0.27,0.73). Inset figure is 
residual Hubble diagram and models after subtracting model 
(f2 m ,,f2A)=(0.27,0.73). The pink line is the luminosity distance 
of the dynamical dark energy equation of state model with (too, 
ui a )=(-2, 8) discussed in § 3-3 



are the luminosity distances of ACDM model with (fl m> 
S1a)=(0.27, 0.73) and dynamical dark energy model with 
(wo, w a )=(-2, 8), respectively. 



3.1 ACDM model 

We first consider the cosmological constant model, that is 
wo — — 1 and w a — 0. Then Eq. Q becomes 



F(z) 



dz 



o y/tt m (i + z) 3 + n A -n k {i + z) 2 
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and x 2 is a function of (fi m ,fiA)- Fig. [2] shows confidence 
regions for (S7„,,S7a) from 63 GRBs (light blue dash-dotted 
lines), 192 SNela (blue dotted lines), and 63 GRBs + 192 
SNela (red solid lines), respectively. Without any prior the 
set of the cosmological parameters with the largest like- 
lihood is (£7 m ,S7 A ) = (0.381^9, 0.93lo;lg) with 
225.2/253. 
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where fi th (zi, fi m , fix, w) = 5 log(d*/Mpc) + 25 and xLst 
represents the chi-square value for the best fit parameter set 
of fl m , Six and w. 



3 COSMOLOGICAL PARAMETERS 

Fig. Q] shows the Hubble diagram extended to z = 5.6 
by GRBs. The green points and red p oints are the lu- 
minosity distance determ i ned by SNela dRiess et al.1 120071 : 
IWood-Vasev et al] 120071 : iDavis et al.1 I2007T I and GRBs 
1 Kodama et al. 20081 ). respectively. The blue and pink lines 



3.2 non-dynamical dark energy model 

In this section we assume a flat universe with w a — 0. Then 



F(z) = 



dz 



y/Qmil + Z) 3 + (1 - HmKi + z) 3(1+ ' 
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and the likelihood function depends on Q m and wo- Fig. [3] 
shows the likelihood contours on (Sl m ,ii>o) plane for GRBs 
(light blue dash-dotted lines), SNela (blue dotted lines), 
SNela + GRBs (solid red lines), respectively. The contours 
correspond to 68.3% and 99.7% confidence regions, respec- 
tively. The set of the cosmological parameters with the 
largest likelihood is (Sl m , wo) — (0.361°'^, —1.33' 
XL* = 227.0/253. 
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Figure 2. The contours of likelihood A% 2 in (fl m fl^) plane 
for GRBs (light blue dash-doted lines), SNeIa(blue dotted lines), 
SNela + GRBs (red solid lines), respectively. The contours cor- 
respond to 68.3% and 99.7% confidence regions, respectively, and 
black solid line represents the fiat universe. The shape of GRB 
contour is more vertical because GRBs have higher redshifts. 
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Figure 4. The contours of likelihood of A\ 2 in (wQ,w a ) plane 
for GRBs (light blue dash-doted lines), SNela (blue dotted lines), 
SNela + GRBs (red solid lines), respectively. The contours cor- 
respond to 68.3% and 99.7 % confidence regions, respectively. 
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Figure 3. The contours of likelihood of Ax 2 in (Q m , wo) plane 
for GRBs (light blue dash-dotted lines), SNela (blue dotted lines), 
SNela + GRBs (red solid lines), respectively. The contours cor- 
respond to 68.3% and 99.7% confidence regions, respectively. 



3.3 Dynamical dark energy model 

As already shown, we adopt the paramete rization of w(z) as 
l|Chevallier fc Polarski|[200ll ; lLinderll2003D 



w(z) = w + w a (l — a) = Wo + w a 



1 + 2 



(7) 



This is not the only parameterization and, for example, 
w(z) — wo + Wiz can be found in the literature. However, 
w(z) is diverging for large z in this (iOo,Wi) parameteriza- 
tion, which would not be appropriate for our high redshift 
GRB samples. Now Eq. © becomes 



F(z) = 



dz 



n m (i + zf 



(l-OnXl+z) 



3 ( 1 + W q + W a ) — 3w a 



-1/2 
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For simplicity, we fix f2 m = 0.28. Fig.|4]shows the contours of 
likelihood A^ 2 in the (too, w a ) plane from GRBs (light blue 
dash-dotted lines), SNela (blue dotted lines), SNela + GRBs 
(red solid lines), respectively. The contours correspond to 
68.3% and 99.7% confidence regions, respectively. The shape 
of probability contour is more horizontal than that of SNela, 
because of the higher redshift distribution of GRBs. The set 
of the dark-energy parameters with the largest likelihood 



is (w , w a ) = (- 



1 26+ 038 1 4+ 10 



with xLst = 227.6/253. 



The pink line in the inset of Fig.[T]is the model with (wo, w a ) 
= (—2,8) which is three sigma level from the best fit. We 
can see that this model does not fit the data by eyes. 

Figs. [2] [3] and [4] show that the constraints on cosmo- 
logical parameters from SNela are stronger than those from 
GRBs at present. However the shapes of the likelihood con- 
tours are different. In Figs.[2]and[3] the contours from GRBs 
are more vertical than that from SNela. This behaviour is 
clear from Eqs. (JS| and (0. GRBs include higher redshifts 
up to z = 5.6 at present so that the value of F(z) is more 
sensitive on the value of Q m . Therefore we can expect inde- 
pendent stronger constraints on cosmological parameters if 
the systematic error in Yonetoku relation decreases and/or 
the number of high redshift GRBs increases. Fig. U shows 
also that the contour from GRBs is more horizontal than 
SNela. Since the mean value of the redshift for SNela sam- 
ples and GRB samples are 0.48 and 1.97, respectively, we 
see that GRBs should give more stronger constraints on w a 
from the exponential term in Eq. §8§ . 



4 FUTURE PROSPECT OF GAMMA-RAY 
BURST COSMOLOGY 

In this section we investigate the future prospect of probing 
dark-energy parameters with GRBs. The Gamma-ray Large 
Area Space Telescope (GLAST) was launched June 11, 2008 
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and it would substantially increase the potential of GRBs as 
cosmological probes. In fact, due to the wide energy-band of 
GLAST Burst Monitor (GBM) and positional accuracy of 
Large Area Telescope (LAT), GLAST is expected to detect 
30 GRBs/year with spectral peak energy and spectroscopic 
redshift by joint observation with Swift . 

Here we estimate the accuracy of determination of the 
dark-energy parameters with GLAST by Monte-Carlo sim- 
ulation. We generate 150 GRB events in the following way. 
First, GRBs are distributed in redshift-luminosity plane ac- 
cording to the GRB formati o n rat e and luminosity function 
given bv lPorciani fc Madaul |200ll) . Spectral peak energy is 
assigned to each GRB according to the best-fit Yonetoku re- 
lation with intrinsic dispersion of 30%. Calculating the flux 
at the earth assuming the concordance model, a GRB is 
counted as an observed event if the flux exceeds the sensi- 
tivity of GLAST. For observed events, observational errors 
of 20% are added to the observed flux and spectral peak 
energy. In this way, we generate 150 observed events and 
they are divided into two groups, high-redshift group (77 
GRBs with z > 1.7) and low-redshift group (73 GRBs with 
z < 1.7). As we did with real GRB events, we reconstruct 
Yonetoku relation with low-redshift GRBs. With increased 
number of low-redshift events, the normalization and index 
of Yonetoku relation are determined with reduced errors of 
5% and 3%, respectively. Applying the reconstructed rela- 
tion to high-redshift events, we can put them in the Hub- 
ble diagram and constrain the cosmological parameters. For 
the details of our Monte-Carlo m ethod, see ((Takahash i et al] 
l2003l ; IOguri fc Takahash"lll2006l ). 

Figs. [S] and HJ] show the constraint in (fl m ,wo) and 
(wo,w a ) planes, respectively, from the "real + simulated" 
high-redshift GRBs (light blue dash-dotted lines), "real" 
SNela (blue dotted lines), and SNela + GRBs (red solid 
lines). The contours correspond to 68.3% and 99.7 % confi- 
dence regions, respectively. These figures indicate how GRBs 
can be a powerful probe to study the nature of dark energy. 

In conclusion, the increase of high redshift GRB data 
by such as GLAST and Swift is indispensable to determine 
the time variation of the dark energy in z > 1.8 where SNela 
data would be rare. 
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fidence regions, respectively. 



6 I r 




-2.0 -1.5 -1.0 -0.5 

wo 
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